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Abstract 

Grosswald’s conjecture is that g{p), the least primitive root modulo p, satisfies g(p) < 
y/p — 2 for all p > 409. We make progress towards this conjecture by proving that 
gip) ^ y/p — 2 for all 409 < p < 2.5 x 10^® and for all p > 3.67 x 10^^. 
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1 Introduction 

Let g{p) denote the least primitive root of a prime p. Burgess [4] showed that g{p) <C 
for any e > 0. This remains the best known bound in general — see [12] for an insightful 
survey of related problems. Grosswald [7] conjectured that 

5'(p)<a/p-2, (1) 

for all primes p > 409. This has implications for the generators of r(p), the principal congru¬ 
ence subgroup modulo p of the modular group T — see [7, §8]. Grosswald verified numerically 
that (1) is true for all 409 < p < 10000. He also gave an explicit version of Burgess’ bound, 
thereby proving that g{p) < pO.499 £qp all p > 1 -|- exp(exp(24)) ~ 10^^ . 

Using computational and theoretical arguments we improve on Grosswald’s estimate in 
the following theorem. 

Theorem 1. Let g{p) denote the least primitive root modulo p. Then g{p) < y/p — 2 for all 
409 < p < 2.5 X 10^® and for all p > 3.67 x lO^U 
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The ‘gap’ in Theorem 1 between the ranges of p seems difficult to bridge. The trivial 
bound g{p) < p when combined with the results in Theorem 1 gives the following corollary. 

Corollary 1. g{p) < for all p. 

The bound in Corollary 1, while weak, appears to be the first bound that holds for all p. 
The remainder of the paper is organised as follows. In §2 we collect the necessary results 
to make Burgess’ result explicit. This gives a substantial improvement on the upper bound 
exp(exp(24)) + 1 given by Grosswald. We introduce a sieving inequality in §3 which enables 
us to reduce this further. Finally, in §4 we present some computational arguments which 
complete the proof of Theorem 1, and present some data on two related problems involving 
primitive roots. 


2 Explicit versions of Burgess’ bounds 


Burgess’s bounds on the character sum 

N+H 

Sh{N)= Y. x{m) 

m=N+l 

were first made explicit by Grosswald [op. cit.], and were later refined by Booker [3], Mc- 
Gown [10], and, most recently, by Trevino [15]. The following is Theorem 1.7 in [15]. 

Theorem 2. [Trevino] Suppose x is a non-principal Dirichlet character modulo p where 
p > 10^^. Let N,H £ 1, with H > 1. Fix a positive integer r > 2. Then there exists a 
computable constant C{r) such that whenever H < we have 

\SHiN)\ <C{r)H^-^/^p'^{logp)^. ( 2 ) 

We follow Burgess, who, in [4, §6] considers 


fix) 


(fjp - 1 ) 

p-l 


1 + 


E 

d\p—l^d>\ 


hid) 

(pid) 


Xd 



(3) 


whence it follows that fix) = 1 if x is a primitive root, and fix) = 0 otherwise. Thus, if 
N(H) denotes the number of primitive roots in the interval N 1 < x < N H we have 


N-CH 

NiH) = Y fi^y 

x=N-el,x^0 (mod p) 


Hence, if H < p there is at most one x £ [V + 1, V + i7] with x = 0 (mod p) so that 


N-CH 

N{H)- Y fi^) 

x=N-\-\ 


< 1. 


We can estimate the sum of /(x) using (2) with H = (1— ^jP^ ■ This choice of H guarantees 

^ Po ^ 

that H < y/p — 2 for p > pQ. 

Since we need only consider square-free divisors d in the outer sum in (3), and since there 
are 4 >id) characters Xd: we arrive at the following theorem. 


2 



( 4 ) 


Theorem 3. TTe have g{p) < ^ — 2 for p > po provided that 

> C{r)(l - -^) '^"(logp)V2»-{2-(P-i) - 1 }. 

^ Po ^ 

The exponent on the left side of (4) is maximised when r = 2. We rearrange (4) accordingly 
to show that we require 

,, ^ , ■ > C(2)i®(0.99)-® 1- 1 (P > 10^°). (5) 

(logpj^ t > 

With C{2) = 3.5751 as in [15, Table 3], we see that (5) is true whenever ijj{p — 1) > 17984. 
Hence we need only consider oj{p — 1) < 17983. Solving for p in (5) we find we need only 
consider p < which is much less than 10^® . We reduce this upper bound substantially 

by introducing a sieving inequality in the next section. 

3 A sieving inequality 

Let e be an even divisor of p — 1. Let Rad(n) denote the product of the distinct prime divisors 
of n. If Rad(e) = Rad(p — 1), then set s = 0 and <5 = 1. Otherwise, if Rad(e) < Rad(p — 1), 
let pi,... ,ps, s > 1, be the primes dividing p — 1 but not e and set 5 = 1 — Yli=iP7^- 
practice, it is essential to choose e so that <5 > 0. 

Again let e be an even divisor of p — 1. An integer x (indivisible by p) will be called e-free 
if, for any divisor d of e, (with d > 1), the congruence x = (mod p) is insoluble. With 
this terminology, a primitive root is (p — l)-free. Given N and H let Ne{H) be the number 
of integers x in the range N-\-\<x<N-\-H that are indivisible by p and such that x is 
e-free. 

Lemma 1. Suppose e is an even divisor of p — 1. Then, in the above notation, 

S 

Np.,{H) > Y^NpAH) -{s- l)Ne{H). (6) 

i=l 

Hence 

S 

Np-i{H) > Yi^pAH) - eipi)Ne{H)] + SNeiH). (7) 

i=l 

Proof. For a given e-free integer x, the right side of (6) contributes 1 if x is additionally 
Pi-free, and otherwise contributes a non-positive quantity. We deduce (7) by rearranging (6) 
bearing in mind the definitions for 6 {pi) and 6 . □ 

Given the divisor e of p — 1, we extend the definition of f{x) to fe{x)^ where 

/e(x) = 0(e) I 1 + 

d\e,d>l ' Xd ) 

and where 0(e) = Hence fe{x) = 1 if x is e-free, and fdx) = 0 otherwise. Thus, 

N+H 

Ne{H) = Y 

x=N+l,x^0 (mod p) 
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It follows from Theorem 2 that, under the constraints of that theorem, 

Ne{H) > e{e)(H- {W{e) - (logp)^ 


where W (e) = is the number of square-free divisors of e. 

Similarly, for any prime divisor / of p — 1 not dividing e, 

Nie{H) - (l - y) ^eiH) < e{e)W{e)C{r)H^-^/^p'^{logp)^, (9) 

where the factor IT(e) arises from the expression W{le) — IT(e). 

Now apply (8) and (9) to (7) to obtain 

Np-i{H) > 5e{e)H - C{r)H^-^/^p'^ {\ogp)^rW{e) + ) • 

V i=i^V 

Since = s — 1 + <5, this yields 

' ^ t 1 pj^ 

Np-i{H) > 6d(e) |i7 - W{e)C{r)H^-^/^p'^ {logp)^r ^ ^) } ‘ 

As in §2, we take H = (1—and r = 2 in (10). This proves the following refinement 
^ Po ^ 

of Theorem 3. 

Theorem 4. Let e he an even divisor of p — 1 and s,S as in Lemma 1 with <5 > 0. ITe have 
dip) < y/P ~ 2 for p > Po provided that 


Oie) 1 - — 1 - C(2) 1 - 


(logp)^ 


1 / s — 1 


-I- 2 ) W{e)p 16 > > 0. (11) 


We can rearrange (11) to show that our criterion becomes 

We consider (12) for uj{p — 1) = n < 17983. By making the choice of s for n given in Table 1 
we verify (12) for all n > 42. 

We are left with those p satisfying w(p — 1) < 41. When uj{p — 1) = n = 41 we choose 
s = 37 to minimise the right-side of (12). This shows that Grosswald’s conjecture is satisfied 
provided that 

> 4.97 X 10®2. (13) 

log p 

Solving (13) for p gives p > 3.67 x 10^^. It is tempting to try to remove the io{p — 1) = 41 
case by enumerating possible primes as in [5]. Since p — 1 > pi • • • p 4 i we seek the number of 
solutions of 

2.98 X 10^° < p < 3.67 X 10^\ p prime, a;(p-l)=41. (14) 

A quick computer check shows that there are 329 different primes that could appear in the 
factorisation of p — 1. While it may be possible to enumerate all such products satisfying 
(14), this would, at best, eliminate the n = 41 case only. We have not pursued such an 
enumeration. 
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Table 1: Choices of s for various ranges of uj{p — 1) = n such that (12) holds. 


Range of u:{p — 1) = n 

s 

[800,17983] 

750 

[400, 799] 

300 

[200, 399] 

180 

[105,199] 

105 

[72,104] 

68 

[55, 71] 

52 

[47,54] 

44 

[43,46] 

40 

42 

38 


4 Computational results 

The computational part of Theorem 1 was proved in the following way. The interval [2,10^®] 
was subdivided into consecutive sub-intervals of manageable size (each with 2^® integers). 
An efficient segmented Eratosthenes sieve (see [2] and [13, §1.1]) was then used to identify 
all primes in each interval. For each prime p that was found, a second Eratosthenes sieve, 
modified to yield complete factorizations [6, §3.2.4], was used to find the factorization of 
{p — l)/2. Since the least primitive root modulo p cannot be of the form with a > 0 and 
6 > 1, i.e., it cannot be a perfect power, the integers 2,3, 5,6, 7,10,... , were tried one at a 
time until a primitive root was found. 

With c as a candidate primitive root, the first test was to check if = — 1 (mod p). 

This was efficiently done using the quadratic reciprocity law data from known tables. If this 
test failed the next c candidate was tried. Otherwise, for each odd prime factor q oi {p— l)/2 
it was checked whether ^ 1 (mod p). The next c candidate was tried if one of these 

tests failed. These tests were efficiently done by performing all modular arithmetic using the 
Montgomery method [11]. Since the “probability” of failure of an individual test is 1/g, the 
odd factors q were sorted in increasing order before performing these tests. Note that g{p) is 
equal to the first c that passes all tests. 

Instead of checking (1) directly for each prime up to 2.5 x 10^®, the record-holder values 
of ^(p), i.e., values of g{p) such that g{p') < g{p) for all p' < p, were computed, as these 
are of independent interest [I] and can be used to check (I) indirectly. The computation 
required a total time of about 3 one-core years, and took about one month to finish on nine 
computers (each with 4 cores) of one computer lab of the Electronics, Telecommunications, 
and Informatics Department of the University of Aveiro. Table 2 presents all g{p) record- 
holders that were found up to 2.5 x 10^^. It extends and corrects one entry of Table 2 of [14], 
which is a summary of computations up to 4 x 10^®. 

The largest g{p) record-holder in Table 2 that does not satisfy (1) is 21, corresponding 
to p = 409. Thus, up to 2.5 x 10^^, the largest p for which (1) is possibly false satisfies 
y/p — 2 < 21, i.e., p < 529. It turns out, as already verified by Grosswald, that the last failure 
of (1) occurs for p = 409. 

An analysis similar to the one described above was also performed for least prime primitive 
roots 5 (p), and for least negative primitive roots h{p). The least negative primitive root 
modulo p is equal to the negative integer, least in absolute value, that is a primitive root 
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Table 2: g{p) record-holders with p < 2.5 x 10^^. 


9{p) 

p 

g{p) 

P 

9{p) 

P 

2 

3 

69 

110881 

179 

6064561441 

3 

7 

73 

760321 

194 

7111268641 

5 

23 

94 

5109721 

197 

9470788801 

6 

41 

97 

17551561 

227 

28725635761 

7 

71 

101 

29418841 

229 

108709927561 

19 

191 

107 

33358081 

263 

386681163961 

21 

409 

111 

45024841 

281 

1990614824641 

23 

2161 

113 

90441961 

293 

44384069747161 

31 

5881 

127 

184254841 

335 

89637484042681 

37 

36721 

137 

324013369 

347 

358973066123281 

38 

55441 

151 

831143041 

359 

2069304073407481 

44 

71761 

164 

1685283601 




modulo p. It cannot be of the form —a^ with a > 0 and b > 2, and is equal to —g{p) if p = 1 
(mod 4). It was found that g{p) < y/p — 2 for 2791 < p < 2.5 x 10^^, and that —h{p) < y/p — 2 
for 409 < p < 10^^. We remark that little is known about either g(p) or h{p) — the reader 
may consult [9] for more details. 


5 Conclusion 

It appears difficult to resolve completely Grosswald’s conjecture. Table 3 in [15] indicates 
that one may hope to reduce the size of C{2) further by taking a larger value of pq. However, 
this appears at present not to give an improvement for our purposes. 

An alternative approach is to use a smoothed version of Burgess’ bounds, in the same way 
that a smoothed Polya-Vinogradov inequality was used in [8]. 
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